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ABSTRACT 



We report on the discovery of a substellar companion or a massive Jupiter orbiting the G5V star HD 16760 with the spectrograph 
SOPHIE installed on the OHP 1.93-m telescope. Characteristics and performances of the spectrograph are presented, as well as the 
SOPHIE exoplanet consortium program. With a minimum mass of 14.3 Mj up , an orbital period of 465 days and an eccentricity of 
0.067, HD 16760b seems to be located just at the end of the mass distribution of giant planets, close to planet/brown-dwarf transition. 
Its quite circular orbit supports a formation in a gaseous protoplanetary disk. 

Key words, planetary systems - Techniques: radial velocities - stars: individual: HD 16760 
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1. Introduction 

The vast majority of 350 known exoplanets have been found 
thanks to radial velocity measurements. Far to be an old- 
fashioned technique, Doppler measurements illustrated these 
last years their capabilities to extend the exoplanet search around 
a large variety of stars. The sensitivity of this technique contin- 
uously increases, opening the possibility to explore the domain 
of low-mass planets down to few Earth masses, to discover and 
characterize multiple planetary systems, to perform long term 
surveys to find true Jupiter-like planets, to establish the plane- 
tary nature and to characterize the transiting candidates of pho- 
tometric surveys. Doppler surveys for exoplanet search require 
high-precision spectrographs and a significant amount of tele- 
scope time over a long duration. 

The SO PHIE spectrograph (Bouchy et al. l2006l Perruchot 
et al. l2008b is in operation since October 2006 at the 1 .93-m tele- 
scope of Observatoire de Haute-Provence. Benefiting from ex- 
perience acquired on HARPS (Pepe et al. 2002) and taking into 
account the limitations of the ELODIE spectrograph (Baranne 
etal.H99(D, SOPHIE was designed to obtain precise radial ve- 
locities with much higher throughput than its predecessor and 
to be operated as a northern counterpart of HARPS. This instru- 
ment is briefly described in section 2. The SOPHIE consortium 
started in October 2006 a large and comprehensive program of 
search and characterization of exoplanets described in section 3. 
We report in section 4 the detection of a substellar companion or 



* Based on observations made with SOPHIE spectrograph on the 
1.93-m telescope at Observatoire de Haute-Provence (CNRS/OAMP), 
France (program 07A.PNP.CONS). 



a massive Jupiter around HD 16760 and we discuss in section 5 
the properties and nature of this object located at the upper limit 
of the mass distribution of giant planets. 

2. The SOPHIE spectrograph 

SOPHIE architecture mainly benefits from ELODIE and 
HARPS experiences. A detailed technical description of this in- 
strument is given by Perruchot et al. (|2008). In this section we 
briefly describe the main properties of the spectrograph and its 
different observing modes. S OPHIE is a cross-dispersed, envi- 
ronmentally stabilized echelle spectrograph dedicated to high- 
precision radial velocity measurements. The detector (EEV- 
4482) is a thinned, back-illuminated, anti-reflection coated 4k 
x 2k 15-/mi-pixel CCD cooled at -100°C, with slow- and fast- 
readout modes. It records 39 spectral orders covering the wave- 
length domain from 3872 to 6943 A. The spectrograph is fed 
through a pair of 3"-wide optical fibers for the high-resolution 
mode (R = 75000, obtained from an extra slit), and another pair 
for the high-efficiency mode (R = 40000, allowing one magni- 
tude gain). The high-resolution mode is equipped with a double 
fiber scrambler (Brown 1990) to homogenize and stabilize the 
illumination of the spectrograph entrance. For each fiber pair, 
one aperture is used for starlight whereas the other one, 2' away 
from the first one, can be used either on a Thorium- Argon lamp 
for tracking spectrograph drift {thosimult mode), or on the sky 
to estimate background pollution, especially in case of strong 
moonlight (objAB mode). Both aperture can also be simultane- 
ously put on Thorium- Argon or tungsten lamps for wavelength 
or flat-field calibrations, respectively. Apart from thermal pre- 
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cautions, the key-point for stability is the encapsulation of the 
dispersive components in a constant pressure tank. This solution 
stabilizes the air refractive index sensitive to atmospheric pres- 
sure variations. With such a concept typical intrinsic drift of the 
spectrograph is less than 3 m s -1 per hour. The ELODIE front- 
end adaptor (Baranne et al. 11996b , is still used for SOPHIE. It 
holds the calibration lamps, the atmospheric dispersion correc- 
tor and the guiding system. Compared to ELODIE, SOPHIE 
leads to 1) gain on photon efficiency by a factor of 10 in high- 
efficiency mode, 2) increase the spectrograph radial velocity sta- 
bility by a factor of 3, and 3) increase spectral resolution from 
42000 to 75000 for the high-resolution mode. 

The spectra are extracted from the detector images and the 
radial velocities are measured online with the S OPHIE pipeline 
derived and adapted from the HARPS on^B The spectra extrac- 
tion includes localization of the 39 spectral orders on the 2D- 
images, optimal order extraction, cosmic -ray rejection, wave- 
length calibration and spectral flat-field correction yielding a two 
dimension spectra (E2DS). The orders are then merged and re- 
binned after correction of the blaze function, yielding a one di- 
mension spectra (SID). The E2DS spectra are cross-correlated 
with numerical masks corresponding to different spectral types 
(F0, G2, K0, K5, M4); the resulting cross-correlation func- 
tions (CCFs) are fitted by Gaussians to get the radial velocities 
(Baranne et al. [19961 Pepe et al. 12002b . 

Following the approach of Santos et al. (2002), we calibrated 
the CCF to determine the projected rotational velocity v sin i and 
the metallicity index [Fe/H] . We also calibrated the CCF to com- 
pute the RV photon-noise uncertainty cr VR . Following the ap- 
proach of Santos et al. (2000), we computed and calibrated the 
chromospheric-activity index R' HK based on our SOPHIE spec- 
tra. 

The SOPHIE radial velocity measurements were initially 
affected by a systematic effect at low signal-to-noise ratio, due 
to CCD charge transfer inefficiency, which increases at low 
flux level. This effect was calibrated and is now corrected 
by the pipeline, which reallocates the charge lost during the 
readout process on each extracted pixel (Bouchy et al. [2009). 
Uncertainties on the radial velocity measurements include pho- 
ton noise, uncertainties in the wavelength calibration, and sys- 
tematic instrumental errors. The photon noise RV uncertainty 
depends on the signal-to-noise of the spectra, as well as on the 
spectral type and the rotation velocity vsinz of the observed 
star. It can be approximated by the semi-empirical estimator 
cr RV = A x VFWHM/(S/N x C), where FWHM is the full width 
at half maximum of the CCF (in same unit as cr RV ), C is its con- 
trast (in percent of the continuum), S/N is the signal-to-noise 
ratio per pixel at 550 nm, and the scaling factor A = 1 .7 or 3.4 in 
high-resolution or high-efficiency mode, respectively. For a non- 
rotating K-dwarf star, a S/N per pixel of 150 provides a photon- 
noise RV uncertainty of 1 m s _1 . Such a S/N is obtained on 5-mn 
on a 6.5 magnitude star. The uncertainty of the wavelength cali- 
bration was estimated to 1 ms _1 . Telescope guiding and center- 
ing errors in averaged weather conditions are typically of 0.3 - 1 
arcsec. In high-resolution mode, these errors imply a RV jitter of 
3-4 ms~' due to the insufficient scrambling gain of the fiber. 
This corresponds to the dispersion obtained on the SOPHIE 
measurements around the orbit of HD 189733b, after correction 
of the stellar jitter (Boisse et al. 2009). Uncertainties due to guid- 
ing errors are more than twice this level in high-efficiency mode 
due to the absence of scrambler in this instrumental setup. 



1 http://www.eso.org/sci/facilities/lasilla/instruments/harps/doc/index.: 



The present radial velocity precision obtained on stable stars 
is about 4-5 ms _1 over several semesters. This limitation is 
mainly due to guiding and centering effects on the fiber entrance 
at the telescope focal plan and the insufficient scrambling pro- 
vided by the fiber and the double scrambler. An upgrade of the 
Cassegrain fiber adapter is presently on-going, including new 
high-precision guiding camera and new double scrambler, with 
the goal to reach the precision level of 1-2 m/s. 

3. The SOPHIE exoplanet program 

The SOPHIE consortium program is devoted exclusively to 
the study and characterization of exoplanets, in continuation of 
a planet-search program initiated 15 years ago with ELODIE 
spectrograph (Queloz et al. 1998) and in complement to the 
HARPS program performed in the southern hemisphere (Mayor 
et al. 120031 1. We started on October 2006 a key program with 
the aims to cover a large part of the exoplanetary science and to 
bring constraints on the formation and evolution processes of 
planetary systems. Our observing strategies and target samples 
are optimized to achieve a variety of science goals and to solve 
several important issues: 1) mass function of planets below the 
mass of Saturn, 2) planetary statistical properties to constrain 
the formation and evolution models, 3) relationships between 
planets and the physical and chemical properties of their 
stars, 4) detection of exoplanets around nearby stars, allowing 
space and ground-based follow-up, 5) deep characterization 
of known transiting exoplanets including long term follow-up 
and spectroscopic transit analysis. All these aspects are treated 
through 5 complementary sub-programs discussed below and 
using an amount of 60 to 90 nights per semester allocated on 
S OPHIE at the 1.93-m telescope. 

- High precision search for super-Earth 

Only a few percents of the 350 detected planets have masses 
less than 0.1 Mj up , and due to the present precision of radial 
velocity surveys, the distribution of planetary masses is heavily 
biased against low-mass planets. Recent HARPS discoveries 
indicate that these low-mass exoplanets are not rare and suggest 
that 30% of non-active G and K dwarfs solar-type harbor 
Neptune or rocky planets with periods shorter than 50 days 
(Lovis et al. 120091 Mayor et al. |2DD9l . From the ELODIE survey 
and from our volume-limited sub-program, we pre-selected a 
sample of about 200 non-active bright solar-type stars to explore 
this domain of low-mass planets. 

- Giant planets survey on a volume-limited sample 

For a large volume-limited sample of 2000 stars, we perform 
a first screening to identify new Hot Jupiters and other Jovian- 
type planets orbiting near and bright stars. Increasing the list of 
Hot Jupiters offers a chance to find a transiting one orbiting a 
bright star appropriate for additional study of planetary atmo- 
sphere. This survey will also provide better statistics to search 
for new properties of the distribution of exoplanet parameters. 
We include on this sub-program the long term follow-up and 
the spin-orbit analysis - from the Rossiter-McLaughlin effect 

- of known transiting giant exoplanets to respectively detect 
additional companions and determine the spin-orbit angle of the 
system. 

- Search for exoplanets around M-dwarfs 

A systematic search for planets is made for a volume-limited 
sample of 180 M-dwarfs closer than 12 parsecs. Such a survey 
«af low mass stars will give us a chance to derive the frequency 
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of planets as function of the stellar mass. The objectives are 
1) to detect exoplanets of few Earth masses in the habitable 
zone, 2) to determine the statistics of planetary systems orbiting 
M-dwarfs in combining these 180-M dwarfs sample with 100-M 
dwarfs monitored with HARPS, 3) to identify new potential 
transiting Hot Neptunes. 

- Search for exoplanets around early-type main sequence stars 

A systematic search for planets around a sample of 300 
early-type main sequence stars (A and F stars) is performed 
to study the impact of the host star mass on the exoplanet 
formation processes. Such stars were previously not included 
in the exoplanet surveys due to their lack of spectral lines and 
high rotation broadening. A specific pipe-line was developed 
to compute radial velocity on these specific targets (Galland et 
al. |2005 a) with an accuracy allowing the detection of planets 
from massive hot Jupiters for fast rotating A stars, and down to 
Neptune-mass planets for the slowest F stars. 

- Long term follow-up of ELODIE long period candidates 

The ELODIE program for exoplanet search which started on 
1994 was performed on a sample of 320 G and K stars. About 
40 of these stars present evidence of long term trends which 
may be due to giant planets with Jupiter or Saturn like orbit. A 
long term follow-up of these candidates is performed to explore 
the domain of long period (> 10 years) planets. 

As part of the S OPHIE consortium programs, the detection 
of four exoplanets have been published up to now: HD 43691b 
an d HP 132406b (Da Silva et al. |2008jl, H P 45652b (Santos et 
al. l2008l . and 9 Cygnib (Pesort et al.|2009). These planets have 
respectively minimum masses of 2.5, 5.6, 0.5 and 2.3 Mj up with 
a 37, 975, 44 and 154 day periods. There were first found from 
the ELOPIE or CORALIE survey then monitored by S OPHIE. 
Spectroscopic transits of the massive planets HP 147506b and 
XO-3b were also observed (Loeillet et al. 2008 and Hebrard et al. 
2008 respectively), allowing a refinement of the parameters of 
the systems, and the detection of a first case of misaligned spin- 
orbit for XO-3 (recently confirmed by Winn et al. 2009). A study 
of the stellar activity of the transiting planet host star HP 189733 
is also presented by Boisse et al. ( I2009I I. Recently the transit 
of the 1 1 1-day period exoplanet HP 80606b was established by 
Moutou et al. ( |20Q9l >. 

Outside of the consortium programs, S OPHIE plays an ef- 
ficient role in the Poppler follow-up of photometric surveys for 
planetary transits search. It allowed the planetary nature to be 
established for transiting candidates found by SuperWASP (e.g. 
Collie r Cameron et al. 120071 Pollac co et al. I20U81 Hebb et al. 
120081 . by HAT (Bako s et al . [20071 and b y CoR oT space mis- 
sion (e.g. Barge et al. 120081 Bouchy et al. 120081 Moutou et al. 
120081 Peleuil et al. 120081 Rauer et al. |2009l , as well as the pa- 
rameters of these new planets to be characterized, including the 
measurement of the masses. 

In the next section we present the detection of the substel- 
lar companion orbiting HP 16760 as part of our sub-program 2 
"Giant planets survey on a volume-limited sample". 

4. The substellar companion of HD 16760 

4. 7. Stellar properties of HD 16760 

HP 16760 (HIP 12638, BP+37 604) is a G5V star located 50 pc 
away according the Hipparcos parallax. Table Q] summarizes the 
stellar parameters. From spectral analysis of the SOPHIE data 



using the method presented in Santos et al. (2004), we derived 
T eS = 5608 ± 20 K, logg = 4.51 ± 0.10, [Fe/H] = -0.02 ± 0.03, 
and M* = 0.86 ± 0.06 M Q , which agrees with values from lit- 
erature. For the temperature and the mass, the values we adopt 
in Table Q] are compromise values between ours and those ob- 
tained by Nordstrom et al. (120041) (T eS = 5636 K and M t = 
O.91+J}^M ). We derive vsin/ = 2.8 ± 1.0 kms -1 from the 
parameters of the CCF using a calibration similar to those pre- 
sented by Santos et al. (2002), in agreement with the value 
vsin/ = 3 kms -1 from Nordstrom et al. (|2004). The CCF also 
allows the value [Fe/H] = 0.06 ± 0.05 to be measured, in agree- 
ment but less accurate than the metallicity obtained above from 
spectral analysis. This target has a quiet chromosphere with no 
emissions in the Can lines (log/?^ K = -5.0 + 0.1), making it a 
favorable target for planets search from radial velocity measure- 
ments. 

Table 1. Adopted stellar parameters for HP 16760. 



Parameters 


Values 


nty 


8.744 


Spectral type 


G5V 


B- V 


0.71 ± 


0.02 


Distance [pc] 


50 = 


= 7 


pmRA [mas/yr] 


82.8 = 


= 3.1 


pmDEC [mas/yr] 


-110.6 


±3.1 


vsin; [kms -1 ] 


2.8 ± 


1.0 




-5.0: 


hO.l 


[Fe/H] 


-0.02 : 


t0.03 


Ttf [K] 


5620 


+ 30 


logs [cgi] 


4.51 = 


= 0.1 


Mass [M Q ] 


0.88 ± 


0.08 



References 
Nordstrom et al. (2004J 
Hipparcos catalog 
Hipparcos catalog 
Hipparcos catalog 
Hipparcos catalog 
Hipparcos catalog 
this work 
this work 
this work 
Nordstrom et al. ( 20041 & this work 

this work 
Nordstrom et al. ( 20041 & this work 



HP 16760 has a stellar companion, HIP 12635 (Apt [19881 
Sinachopoulos 2007), located 14.562 ± 0.008 arcsec in the 
North and 1.521 ± 0.002 mag fainter. From Hipparcos cata- 
log (Perryman et al., 1997), HIP 12635 has similar distance 
(45 + 17pc) and similar proper motions (pmRA=107 ±17, 
pmPEC=102 ± 12) with HP 16760, making them a likely phys- 
ical system, with a separation > 700 AU and an orbital period 
> 10 000 years. This would induce tiny radial velocity variations 
on the stars, below 0.2ms _1 yr _1 . 

4.2. Radial velocity measurements and Keplerian fit 

We acquired with SOPHIE 20 spectra of HP 1 6760 within objAB 
mode between Pecember 2006 and October 2008 under good 
weather conditions. Two of these 20 spectra were polluted by 
significant Moon contamination. The velocity of the CCF due to 
the Moon was far enough from those of the target to avoid any 
significant effect on the radial velocity measurement. We did not 
use these two spectra however for the spectral analysis presented 
in§E] 

The derived radial velocities are reported in Table[2] together 
with the journal of the observations. Their typical accuracy is 
around 6 ms -1 , which is the quadratic sum of three sources of 
noise: photon noise (3 ms -1 ), guiding (4 ms 4 ), and spectro- 
graph drift (3 ms -1 ). 

The radial velocities, shown in Fig. Q] present clear varia- 
tions of the order of hundreds ms -1 , without significant vari- 
ations (<x < 10 ms -1 ) of the CCF bisector (Fig. thus in 
agreement with the reflex motion due to a companion. We fit- 
ted the data with a Keplerian model. The solution is a 465-day 
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Fig. 1. Top: Radial velocity SOPHIE measurements of 
HD 16760 as a function of time, and Keplerian fit to the data. 
The orbital parameters corresponding to this fit are reported in 
Table[3] Bottom: Residuals of the fit with 1-cr error bars. 



Table 2. Radial velocities of HD 16760 measured with SOPHIE. 



BJD 


RV 


±lo- 


exp. time 


S/N p. pix. 


-2400000 


(kms" 1 ) 


(km s" 1 ) 


(sec) 


(at 550 nm) 


54099.3679 


-3.6736 


0.0055 


300 


62.5 


54126.3439 


-3.7916 


0.0058 


300 


52.1 


54127.2914 


-3.8085 


0.0068 


480 


36.2 


54133.2541 


-3.8478 


0.0055 


224 


61.4 


54138.3414 


-3.8366 


0.0054 


673 


67.1 


54148.3393 


-3.8898 


0.0071 


224 


35.1 


54151.3375 


-3.9150 


0.0092 


225 


30.0 


54155.2737 


-3.9168 


0.0063 


180 


43.5 


54339.6353 


-3.3666 


0.0053 


512 


78.7 


54352.6607 


-3.2786 


0.0058 


380 


52.0 


54367.6375 


-3.2481 


0.0055 


300 


61.2 


54407.4733 


-3.1806 


0.0053 


500 


79.5 


54434.4441 


-3.1703 


0.0055 


620 


64.2 


54502.2489 


-3.3798 


0.0057 


443 


54.5 


54513.2880 


-3.4358 


0.0059 


464 


53.4 


54546.2900 


-3.5725 


0.0061 


1304 


47.3 


54683.6329 


-3.9680 


0.0058 


393 


51.4 


54705.6646 


-3.9045 


0.0058 


379 


51.2 


54739.5899 


-3.7242 


0.0059 


350 


50.1 


54755.5295 


-3.6488 


0.0058 


270 


51.2 



period oscillation with a semi-amplitude K = 408 ms -1 , cor- 
responding to a substellar companion, with a minimum mass 
nip sin i - 14.3 Mj up . The derived orbital parameters are reported 
in Table [3] together with error bars, which were computed from 
X 2 variations and Monte Carlo experiments. 

The standard deviation of the residuals to the fit is cr(O-C) = 
10.1 m s _1 . This is higher than the 6- m s _1 estimated uncertainty 
on the individual measurements. Although about 23% of gaseous 
giant planets are in a multiple planetary system, we do not iden- 
tify yet an indication for a second body orbiting HD 16760. With 
a maximum semi-amplitude of 20 ms _1 , the residuals of the fit 
do not exhibit structures, denying a possible inner planet with a 
rap sin i > 0.7 M j up . A longer period planet may induce a drift 
lower than 20 m s~'yr~' during our observational period. 



0.05 - 



tn 



E 0.00 - 




-0.10 - 



-4.0 -3.8 -3.6 -3.4 -3.2 

Radiol velocity (km/s) 

Fig. 2. Bisector span as a function of the radial velocity. 



Table 3. Fitted orbit and planetary parameters for HD 16760b. 



Parameters 


Values and l-cr error bars 


Unit 


V r 


-3.561 ±0.004 


kms"' 


P 


465.1 ±2.3 


days 


e 


0.067 ±0.010 




CO 


-128 ± 10 




K 


408 ±7 


ms" 1 


T (periastron) 


2454723 ± 12 


BJD 


cr(0 - C) 


10.1 


ms" 1 


reduced \ 2 


2.0 




Nobs 


20 




nip sin i 


14.3 ±0.9" 


M Jup 


a 


1.13±0.03 f 


AU 



t: using M* = 0.88 ± 0.08 M, 



5. Discussion and Conclusion 

Our RV measurements indicates that a substellar companion 
with aprojected mass rap sin i = 14.3 Mj up is orbiting HD 16760. 
With the degeneracy of inclination angle i, it is difficult to con- 
clude about the exact nature of this companion. It may corre- 
spond to a massive planet, formed in a gaseous protoplanetary 
disk, or a brown dwarf, issued from collapse in a giant molecu- 
lar cloud. 

Figure [3] shows the mass distribution of massive planets 
(rap sin i > 3 Mj up ) and light brown dwarfs (Mc sin i < 30 Mj up ) 
found by radial velocity surveys. From the Extrasolar Planets 
Encyclopaedia lisfl, we completed with HD 137510b (Endl et 
al. 120041 . HD 180777b (Galland et al. l2005bt . and HD 16760b 
(this paper), totalizing 89 objects including 10 with mass in- 
between 15 and 30 Mj up . The dashed curve corresponds to the 
relation M~ 2 (dN/dM=M _1 ). The black shaded histogram corre- 
sponds to the transiting planets with true masses (excluding non- 
confirming objects SWEEPS- 11 and SWEEPS-04). It is worth- 
while to notice that, although based on a small number of object, 
the mass distribution of transiting planets is following the same 
trend than non-transiting planets. Indeed, the ratio of transiting 
planets over non-transiting planets is about the same : 9.2, 9.5 
and 10 for the bins 3-6, 6-9 and 9-12 Mj up respectively. In this 
histogram, HD 16760b seems to be located just at the end of the 
mass distribution of giant planets. Although based on small num- 
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bers, sub-stellar companions with minimum mass greater than 
17 Mj up do not seem to follow the M~ 2 relation. 

Figure |4] shows the eccentricity - period diagram of mas- 
sive exoplanets and light brown dwarfs. The size of circle is 
function of the mass (3-5, 5-10, 10-15 Mj up ). Hexagonal points 
corresponds to objets with mass greater than 15 Mj up . Black 
filled symbols correspond to transiting companions. HD 16760b 
confirms the observed trend that more massive companions are 
found for longer period planets (Udry & Santos 2007 ). We also 
notice that all companions with mass greater than 15 Mj up have 
an eccentricity greater than 0.2 except CoRoT-exo-3b (Deleuil et 
al. 120081 and HD41004Bb (Zucker et al. [2004b in very close-in 
orbit (with periods of respectively 4.2 et 1.3 days) tidally cir- 
cularized. The properties of HD 16760b make it an interresting 
sub-stellar companion. With a mass greater than the Deuterium 
burning limit (13 Mj up ), it may be defined as a brow-dwarf. 
However its quite circular orbit supports a formation in a gaseous 
protoplanetary disk. On another way, Halbwachs et al. (2005) 
studied the eccentricity distribution for exoplanets and binary 
stars with a mass ratio smaller than 0.8 (non twin binaries). They 
found that exoplanets have orbits with eccentricities significantly 
smaller than those of the non-twin binaries, reinforcing the hy- 
pothesis that planetary systems and stellar binaries are not the 
products of the same physical process. 

HD 16760b is in a visual and physical binary system. 
However, in the mass-period and eccentricity-period diagrams, 
HD 16760 b is located in a region not much populated by plan- 
ets in binary systems. The discovery of this long-period low- 
eccentricity planet thus adds to the growing evidence that con- 
trary to short-period planets, long-period (>100 days) planets 
residing in binaries possess the same statistical properties as 
their counterparts orbiting single stars (Eggenberger et al. 2004 
Mugrauer et al. 120051 Desidera et al. |20071 >. 

HD 16760b would induce a motion of its host star of at least 
±0.35 milli-arcsec. The future Gaia ESA space mission sched- 
uled for launch in late-2011, should be able to detect this sys- 
tem from astrometry, and thus would allow the inclination of 
the system to be measured and the true mass to be determined. 
Detailed caracterization of this sub-stellar companion close to 
planet/brow-dwarf transition will help to distinguish the differ- 
ences of formation processes between these two populations. 
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